Geometric ergodicity of a bead- spring pair with 
stochastic Stokes forcing 



Jonathan C. Mattingly, Scott A. McKinley, Natesh S. Pillai 
February 25th, 2009 

Abstract 

We consider a simple model for the fluctuating hydrodynamics of a flexi- 
ble polymer in dilute solution, demonstrating geometric ergodicity for a pair 
of particles that interact with each other through a nonlinear spring poten- 
tial while being advected by a stochastic Stokes fluid velocity field. This is 
a generalization of previous models which have used linear spring forces as 
well as white-in-time fluid velocity fields. 

We follow previous work combining control theoretic arguments, Lya- 
punov functions, and hypo-elliptic diffusion theory to prove exponential con- 
vergence via a Harris chain argument. To this, we add the possibility of 
excluding certain "bad" sets in phase space in which the assumptions are 
violated but from which the systems leaves with a controllable probability. 
This allows for the treatment of singular drifts, such as those derived from 
the Lennard- Jones potential, which is an novel feature of this work. 



1 Introduction 

The study of polymer stretching in random fluids has been identified as a first 
step in the much larger project of modeling and understanding drag reduction in 
polymer solutions [CheOd] and theoretical focus has been brought on the dynamics 
of simple dumbbeh models |LMV02l, P CMV05L IIAV05L Of particular interest 
is the experimentally observed phenomenon called the coiled state / stretched state 
phase transition IIGCS05II . Mathematically this transition has been characterized by 
seeking models can admit solutions that are ergodic for certain regions of parameter 
space, while being null recurrent for other parameters |CMV05|. In this paper 
we address the topic of how to prove ergodicity for a wide range of models that 
generalize preceding work. 

Let xi(i) and X2(t) denote the respective positions in of two polymer 



"beads" comiected by a "spring" at time t. As noted in IIDE861 10961 these beads are 
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not intended to model individual monomers, nor is the spring intended to capture 
the mechanics of an actual molecule. Rather the hope is to study qualitative fea- 
tures of generic chains interacting in a randomly fluctuating incompressible fluid 
environment. 

Having made this caveat, the canonical Langevin model for two spherical par- 
ticles in a passive polymer system is given by 

mxi = -V$(|xi -xjl) +C(u(xi(t),t) -Xi(t)) + KW(t) . (1) 

The mass m is considered to be vanishingly small and so the inertial term, mxj, 
will be ignored. On the right hand side, the first term is the restorative force exerted 
on the beads due to the potential energy of the polymer's current configuration. The 
second term is an expression for the drag force exerted by a time-dependent fluid 
velocity field u with friction coefficient := Qirarj following from Stokes drag 
law for a spherical particle of radius a in a fluid with viscosity rj. The final term 
is the force due to thermal fluctuations in the fluid. The diffusive constant k is 
often taken to be k = v^2fc^TC, where A;^ is the Boltzmann constant and T is the 
temperature of the system in Kelvin, in accordance with the fluctuation-dissipation 
theorem IICMV05L 

The goal of the present work is to achieve rigorous results about the ergodicity 
of the connector process 

r(t) := i(xi(t)-X2(t)) 

in both the k = and k 7^ regimes with nonlinear spring interaction in the pres- 
ence of a spatially and temporally correlated incompressible fluid velocity field. 

In the simplest possible setting one ignores the fluid and assumes a Hookean 
(quadratic) spring potential <I>. In this case ([T|l is a simplification of the classical 
Rouse model l,DE86.l . We define 

^(x,y) := -^Ix-yP 

and the particle dynamics satisfy the system of SDE 

dxi (t ) = 7 [x2 ) - XI (t )] + K dWi (t ) 

dx2 [t) = 7 [Xi (t ) - X2 (t )] + K dW2 {t) 

where Wi and W2 are independent standard Brownian motions. The dynamics of 
the connector r(t) are given by 

dT{t) = -7r(t) + K\/2(fW(t). 
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where W = ^(Wi - W2) is a standard Brownian motion. We see that each of 
the connector components is an Ornstein-Uhlenbeck process which therefore has 
the unique invariant measure 



This exactly solvable model does not yield physical results, so one must adopt 
nonlinear models for either or both of the spring potential and fluid forces. 

Significant theoretical advances exist for the dynamics of a single tracer par- 
ticle convected by a wide variety of fluid models ltMK99ll . One popular fluid 
model for non-interacting two-point motions |BCH07] iMWD"'"05l as well as for 
Hookean bead-spring systems llCheOOi ILM V02i, iCMV05]| is a time-dependent ran- 
dom field satisfying the statistics of the Kraichnan-Batchelor ensemble I.Bat59l 
IIKra68i . Such a fluid is still white in time, but is coloured in space. 

In the case where k = with non-interacting beads, the spatial correlations 
in the convecting fluid velocity field allow for concentration and aggregation phe- 
nomena [,SS02b l IMWD+05.1 I BCH07 I . This happens because when the two beads 
are very close together, the fluid forces on the respective beads are so strongly cor- 
related there is no force encouraging separation. 

The presence of a diffusive term with k / prevents such aggregation and 
the long term behavior of the connector depends on so-called Weissenberg number 
Wi = IICMV05II . It is shown that when Wi < 1 the connector r will have 
a non-trivial stationary distribution, dubbed the "coiled" state. For Wi > 1, the 
connector does not have a stationary distribution and is called "stretched." The 
authors express interest in the case where the fluid is not assumed to be white-in- 
time. 

In this work we use the incompressible stochastic Stokes equations to gener- 
ate a fluid that is coloured in space and time (see Section 11.21 ). In the Hookean 
spring case with k = 0, this model leads to degenerate dynamics (Proposition 1 1.1b . 
However, in a more general setting with a non-linear spring potential, we show 
that dynamics are nondegenerate, although the coiled / stretched state dichotomy 
discussed in IICMV05II is not present. We find that r{t) is ergodic regardless of the 
physical parameters (Theorem 12. lb . 

The method used here to establish ergodicity builds on the Harris Chain theory 
developed in I Har561 iHasSOl lNum84i It is particularly indebted to the uniform 
ergodic results in weighted norms developed in [ MT 93al lMT93bll . The argument 
follows the path outlined in IIMS02II IIMSH02II for unique ergodicity of degenerate 
diffusions, but requires some nontrivial extensions to deal with the multiplicative 
nature of the noise and to permit the type of singular vector fields that arise as 
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natural choices for the spring potential We build a framework around a general 
ergodic result from IIHM08II and then develop the needed analysis to apply this 
framework. 

Mathematically, as in IIMSH021 IMS02 1. this paper combines control theory 
with techniques from the theory of hypoelliptic diffusions to invoke results in the 
spirit of ||MT93a[|MT93bll . where ergodicity is obtained by proving a minorization 
condition on a class of "small sets" (see IIMT93al IMT93b1 ) while establishing a 
matching Lyapunov function. However, the current setting presents a number of 
difficulties which prevent the application of the results |MSH02 1 directly. In par- 
ticular the spring potential, and hence the drift term, has a singularity (Assumption 
|2ll. Therefore the natural candidates for "small sets" are not compact. This diffi- 
culty is overcome by splitting the small sets into "good" and "bad" sets. On the 
compact "good" set, defined in Eq. ([TT] ). we demonstrate uniform controllability 
as in IIMSH021 IMS02II . On the bad set, one cannot obtain uniform control; how- 
ever, the deterministic dynamics move the system into the good set in finite time so 
that geometric ergodicity still holds (Section [Z2l ). Allowing the spring potential to 
be singular extends the applicability of the theory to many interesting, physically 
important potentials such as the Lennard- Jones potential. Related ideas have been 
also recently been used to prove related, but different, ergodic and homogenization 
results in different settings (see iBub09[|HP07il ). 

1.1 Structure of Paper and overview of results 

For the remainder of Section 1 , we propose the model and explore its dynamics 
when the distance between the two beads, r, is close to zero. Proposition 11.11 
shows that when the spring is Hookean with zero rest length (a quadratic potential), 
r(t) ^ as f ^ oo almost surely if the spring constant 7 is sufficiently strong 
relative to a quantity that depends on the typical spatial gradients in the random 
forcing. Lemma [L2] shows that, in the case of a nonlinear spring with non-zero 
rest length, r{t) never equals zero nor converges to it as t ^ 00. 

In Section [2l we quote an abstract ergodic result from the literature through 
which the results in this paper proceed. The quoted result requires proving a mi- 
norization condition and the existence of a Lyapunov function. Section 12.11 con- 
tains a general prescription for how to deduce the minorization condition from 
the existence of a continuous transition density and a weak form of topological 
irreducibility for the Markov process. In Section 12.21 the needed topological ir- 
reducibility is proven via a control theoretic argument. In Section 12.31 we invoke 
Hormander's "sum of squares" theorem to prove that the associated hypoelliptic 
diffusion has a smooth transition density. Section 1241 contains the calculations es- 
tabhshing the existence of a Lyapunov function and Section l23] contains a number 
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of generalizations and implications of the preceding results. The appendix contains 
the derivation of the model used. 



1.2 Definition of the model 

In the overdamped, highly viscous regime, it is reasonable to neglect the nonlin- 
ear term in Navier-Stokes equations. Following [ OR89.I . liMS02J . IIMSH02 1 and 
IISS02al we consider the bead-spring system advected by a random field u satis- 
fying the incompressible time-dependent stochastic Stokes equations. Following 
IIWal86L IIDZ92L IIDal99ll and IIMcK06ll we have the stochastic PDF 

dtu{x, t) - z^Au(x, t) + Vp(x, t) = F{dx, dt) (2) 
V- u(x,t) = 

with periodic boundary conditions on the rectangle R := [0,L] x [0,L] where L 
is presumed to be very large. We assume that the space-time forcing is a Gaussian 
process with covariance 

E[F(x,t)]=0, K[F\x,t)F^ {y, s)] = {t A s)2kBTiy5,jT i^-y) (3) 

where t A s is the minimum of t and s; the components i, j G {xi,X2} and 6ij is 
a Kronecker delta function; and the remaining constants each have physical mean- 
ing: A;^ is Boltzmann's constant, T is the temperature of the system and v is the 
viscosity. As is shown in the Appendix, it follows that 



F(x,i) = ^^^ Y: e'-^^-/^aMt) . (4) 
kez2\o 

where the are independent standard 2-d Brownian motions and the coefficients 
Uk are related to the spatial correlation function T through the Fourier relation 

fcez2\{o} 

This relation is possible because F is a covariance function, and therefore positive 
definite. By Bochner's Theorem, F is realizable as the the Fourier inverse transform 
of a positive real "spectral" measure. In the periodic setting, this is the counting 
measure 

f (k) = ^j2kBTval. 

Often, one defines the correlation structure on the Fourier side directly. We choose 
a radial, summable shape structure for the {cTk}. 
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Assumption 1. Let ^ : — be given with = <i){\'x\) and let the set 

active mode directions K. G TP be given. We define := ^^(k) and require that 
the norm 

is finite for all s > 1. 

Our interest will be rigorous analysis of the long-term behavior of the connector 
process r whose dynamics, after shght simplification, are given by the following 
system (see Appendix for details). We define the Markov process X(t) = {^^^J^} to 
be the unique solution to the system of SDE 




where the satisfy 

dzk{t) = -AV|kpZk(i)c?t + y/2p^akdWk{t) (6) 

where A = 2ir/L and P = kBT/An^. We define the norm on this family of 
processes, ||z|p := Y^keK l-^kP- 

We note that the eigenmodes {zk} are a mutually independent collection of 
Omstein-Uhlenbeck processes which have the stationary distribution 

and autocorrelation function 

E[zk(t>k(s)]=/?|^e-^'l*-^l . 

For a position r G and the family {^kjkeA: chosen from their respective sta- 
tionary distributions, we note that the state-dependent fluid forcing 

U(r,z) := ^sin(Ak-r)|^Zk (7) 

is a mean zero random variable whose variance is bounded by a constant times the 
Hi norm of the spectral measure 

E[U(r,z)] = 0, E[|U(r,z)|2] < f5\\a\\l < oo. 
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Furthermore, we will use the pathwise estimate 

|U(r,z)|2 < ||zf . (8) 

As mentioned earlier, the choice of quadratic potential $ corresponds to a 
Hookean spring model, but this yields degenerate dynamics (Proposition ll.ll ). In 
the sequel, we place the following assumptions on the nonlinear spring potential. 

Assumption 2. Let $(r) = $(|r|) be a continuously differentiable function satis- 
fying the following 

(i) For every R>0, the set {r G M_(_ s.t. <I>(r) < R} is compact. 

( ii) There exists an Rq such that for all r > Rq 

V$(r) -r > 7|rp (9) 
and there exists c > and eo > such that for all r satisfying |r| G (0, eo] 

- V$(r) • r > c (10) 

Remark 1. In is in the context of this assumption that we choose the length of the 
periodicity of the forcing fluid. Henceforth we take L ^ 4i?o- 

The above description includes as an example the linear spring potential with 
non-zero rest length. In this case, 

^r) = l{\r\-Rf 

for some constant R denoting the rest length. The theory presented also allows for 
potentials with singularities, as seen in the family of functions 

$(r) = -L|r|2'? + ^ (11) 

where q > 1 and a > 0. The choice of g = 1 and a = 12 corresponds to a 
Lennard- Jones singularity at zero and a classical linear spring at infinity. 

A perhaps more common choice is the finite extensible nonlinear elastic FENE 
model 

^(r)=ln(^l-^ 

While this choice is acceptable away from zero, this logarithmic term alone does 
not satisfy the near-zero condition and some repulsive force must be added at the 
origin. 
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1.3 Near-zero dynamics 

Behavior near the origin is critical to determining whether our bead-spring model 
supports non-trivial dynamics. The following quick calculation demonstrates that 
without a repulsive force at the origin, under mild conditions on the spring constant 
and spectral measure of the fluid, the two beads become trapped near each other in 
the long run almost surely. 

Proposition 1.1 (Degeneracy of the Hookean spring case). Let r and the family 
{^kjkeA: satisfy the system of differential equations ^ and (O, with the spectral 
measure {a\^ satisfying Assumption\l\ Let^(r) = ^\r\^ and suppose that \\(j\\o < 
oo. Then there exists a 7o = Todlfllo) tl^t ifl > 7o then 



lim r(t) = 

t~>co 



almost surely. 



Proof. Note that the process |r(i)p satisfies the following pathwise ODE, where 
u) denotes the sample space variable. 



-|r(t; = -27|r(t; a;)|2 + 2 J]] sin(Ak • r(t; a;)) u;) 

ke/c ' ' 

< -27|r(t;..)|2 + 2A J] |k • r(t;u;)||k^ • r(t;^)|^!i|^ 

< -27|r(t; uj)\^ + 2A|r(t; u:)\^ ^ |k| |zk(t; 

kGX: 

This differential inequality implies (and suppressing the dependence on uj) 

ft 



-2jt + 2X [ V |k||zk(s)|ds 
^0 ke^ 



|r(t)|^ < exp 

By the Law of Large Numbers 

^ f Y.MzUs)\ds^y^\\a\\ 



almost surely as t —>■ oo, and so we see that for sufficiently large 7 one has that 
|r(t)p almost surely as i ^ 00. 

□ 
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The remainder of the paper is devoted to the nonlinear spring case. We first 
establish that the origin is an unattainable point for the process r{t) and that after 
a certain interval of time, there is a positive probability of escape from a neigh- 
borhood of the origin. Due to the exponential decay of memory in the system, 
we have in fact shown connector process r will escape a near-origin neighborhood 
with probability 1. 

We begin by noting the following growth inequality. Since the strength of the 
fluid forcing is limited in a neighborhood of the origin, the deterministic spring 
force dominates and sends the two beads apart with positive probability. 

Lemma 1.2. Suppose the spectral measure {cJk} satisfies Assumption \I} and the 
spring potential $ satisfies Assumption |2] Then the origin is unattainable for the 
connector vector process, i.e. |r(t)| > Q for all t almost surely. Furthermore, v{t) 
leaves any sufficiently small open neighborhood of the origin infinite time. 

Proof. Fixing any M > 0, we define to = and tn = inf{t > 1 + tn-i '■ 
||z(t)|| < A/}. Standard properties of the ensure that t„ < oo with probability 
one. Furthermore, for any M > M there exists an q > so that > o 

where 0„ = {supggjj^ t„+i] ^ AT}. Let eo and c be the constants from 

Assumption |2] Then for any -d G (0, 1) and e G (0, eq] if one sets r = inf{s > : 
|r(s)| > e} then for t G [0, r], 

|i|r(t)|2 = -V^.(r(t)) • v{t) + U(r(t),z(t)) • r(t) 

>c-i^|r(t)p-^|U(r(t),z(t))|2 

>c-^\v{t)\'-n^{t)\\. 

where we recall the u-hy-u estimate ([8]). Further restricting to any a; G f^n and 
fixing I? = c/ (2M), one has 

d 1 , , , ,9 c M , , , ,9 
Assuming that r > integrating the preceding estimate produces 

rtAT 

|r(t Ar)|2 > e~(*^^-*")*^/^|r(t„)|2+c / e-^'-'-^^' ds 

J tn 

> ^fl _ g-(lA{r-t„))A^M 

Fixing e = eo A (1 - e-^l'')lM, we see that on J7„ it is impossible to have 
r > tn + 1. Hence we see that P{r G [tn, in + > > Using the strong 
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Markov property of the family z and the fact that the t„ are stopping times, we 
have P{r > tn + 1} < [l — a)" which concludes the proof. □ 

We have in fact proven the following which will be used later. 

Corollary 1.3. Given M > M > 0, there exists an e and an a G (0,1) such that 

r(ro,zo) := inf{t > : |r(t)| > eand ||z(t)|| < M} 

satisfies 

inf inf P{T(ro, Zq) < 1} > a 

{zq: ||zo||<Af} {ro: 0<|ro|<eo} 

2 Ergodicity 

We begin by setting some notation. Let X(t) = (^j*)) satisfy the system Q-®. 
For technical reasons we henceforth restrict ourselves to the case where the set of 
active modes K, (from Assumption [T]) is finite. It follows from Proposition 1 1 .21 that 
the Markov process X(t) is well-defined on the state space 

X:= {(r,z) G (M^ \ Q) x R^} . 

For a bounded, measurable function : X ^ M we define the action of the Markov 
semigroup Vt by 

(Pi^)(x)=Ex[(^(X(t))]. 

To measure convergence to equilibrium we introduce the following weighted 
norm on such functions Lp relative to a given Lyapunov function y : X ^ [0, oo), 

II II l'^^^)! 

:= sup ■ 



1 + V^(x) ' 
In this paper, we take 

y(x) := |r|2 Vi?^ +r/|zp 

where Rq is the value given in equation ^ of the spring potential Assumption |2] 
and ?7 is a constant to be chosen later in Lemma 12.71 We note that the Markov 
semigroup Vt can be extended to act on all functions (p bounded pointwise above 
by V. 

The main result of this article is the following statement about the geometric 
ergodicity of the full Markov process X, and by corollary, the marginal process r 
converges to a unique non-trivial stationary distribution in exponential time. 
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Theorem 2.1. Suppose that the set of active modes /C is finite, but contains at 
least three pairwise linearly independent vectors, and let the spring potential <I> 
satisfy Assumption |2] Then there exists a unique non-trivial invariant measure vr 
and constants C > and \> Q so that 

\\Vtip - TTipW < Ce~^^\\ip\\ 

where irip = J ipdn. 

We begin by recording in a more abstract setting the norm in which the dynam- 
ics converge to equilibrium. It is a weighted version of the total variation distance 
between measures. We generalize the preceding weighted L°°-norm with a family 
of equivalent norms depending on a scale parameter (3 > 0. For any ip: X —>■ W 
define 

We use this to define an associated norm on probability measures. 

P/3(Mi>At2)= sup ip{x.)ni{dx) - ip{x)ij,2{dx) 

The conclusion of Theorem l2. 1 [ follows from the main theorem in Ref IIHM081 
adapted to our setting: 

Theorem 2.2. Suppose that for some t > 0, ci > and cq G (0, 1), a function 
y : X ^ [0, cxd) having compact level sets with lim||x||^oo ^(x) = ^ satisfies 

(7't^)(x) <coF(x) + ci (12) 

for all X G X. Furthermore suppose there exists a probability measure v and 
constant a £ (0, 1) such that 

inf Vt{x, ■ ) > au{ ■ ) (13) 

with C := {x G X: V{x) < R} for some R > 2ci/(l - cq). 
Then there exists an oq £ (0, 1) and (3 > so that 

P/3(^tVl,^t>2) < aoPl3{lll,l^2) 
for any probability measure fii and fi2 on X. 
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The first condition ([T2l ) states that F is a Lyapunov function for the dynamics, 
established in Lemma 12.71 In Lemma 12.61 we construct a minorizing measure, as 
required by condition ([T3] ). This lemma follows from the combination of a form of 
topological irreducibility (Proposition 12.41 ) and local smoothing (Proposition 12.61 ). 
The local smoothing follows from hypoellipticity of the generator of the Markov 
process X and a version of Hormander's sum of squares theorem (cf. IIHor85[ 



2.1 Conditions for measure-theoretic irreducibility 

We now show how to use a very weak topological irreducibility to prove the 
measure-theoretic minorization/irreducibility property given in ([T3l) . We begin by 
fixing the set C which should be thought of as the "center" of the state space: 

C := {x G X : ^(x) < 2i?g}. (14) 

Proposition 2.3. If the following two conditions hold, then there exists a constant 
a G (0, 1) and probability measure v so that (1131 ) holds. 

(i) Uniformly Accessible Neighborhood Condition: There exists a x* G C such 
that for any (5 > there exists a ti > Q and a function ai = ai(t, 6) > 0, 
continuous in t, such that 



mfVt{^,Bs{^,))>ai (15) 
xec 



for all t > ti. 



(ii) Continuous Density Condition: There exists s > and an open set O C C, 
with x^, G O, such that for x G O and measurable A C O one has 



Ps(x, A) = / psix,y)dy 
J A 



with Ps(x, y) jointly continuous in (x, y) for x, y G O and y*) > 

for some G O. 

Proof. By the continuity assumption on ps there exists 5 > so that Bs (x^, ),Bs{y^) C 
Oand 

inf inf Ps(x,y) > ^Ps(x*,y*) > . 

xeB^Cx.) yeB^(y,) Z 

We define the minorizing probability measure u by 

XiAnBsiy,)) 



12 



where A is Lebesgue measure and A is any measurable set. 

For any t > ii + s, we define a{t) = ^ps{x^,y^,)ai{t — s), where ai is the 
function from the Uniformly Accessible Neighborhood Condition Given any 
measurable set A and xq G C we have 

Vtixo,A)= [ [ Vt,ixo,dx.)Vsix,dy) 

> / 'Ptj{xo,dx)ps{x,y)dy > au{A) , 

J AnBsiy,) JBsix,) 

which proves the claim. □ 



2.2 Topological irreducibility via controllability 

Recall the "center" C of the state space X from ([T4l ). In order to demonstrate the 
Uniformly Accessible Neighborhood Condition ^ given in Lemma 12.31 we wish 
to use the z process to drive the r process to some reference point r* . However, due 
to the possible singularity at the origin (see Assumption |2l) the differential equation 
([5]) for r may have unbounded coefficients. We therefore will designate a region of 
bad control within the center C, as well as a compact region of good control. 

To this end, let ei be the constant derived from applying Corollary 11.31 with 
M = i?o and M = V2Ro/ri. We define the set of "bad" points in C by 

^={(r,z)GC:|r|<ei}. (16) 

We define the set of "good" points Q to be the remaining portion of C plus an 
additional region for technical reasons. Precisely, 

g = g^xg^:= |(r,z) G X : |r| G [ei, y2i?o] , ||z|| < m} . (17) 

where the constant M := max{^^, |/Cp(l + maxreg^ |V<I>(r)|)} with r] defined 
later in Lemma IZTl 

We now use a controllability argument to establish the weak form of uniform 
topological irreducibility on g given in Eq. ([T5] ) (for the set C). 

Lemma 2.4 (Topological irreducibility on the "good" set g). There exists a refer- 
ence point X* G C such for any (5 > 0, there exists ati > 50 that for any closed 
interval Id [ti, oo) there is an ai> Q with 

inf inf Pt(x,55(x*)) > ai (18) 
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Sketch of Proof: We begin by finding a suitable reference point x* G Q. By 
Assumption [2j the spring potential $ has a (possibly non-unique) global minimum 
r^,. Since the global minimum of the norm || • || is achieved at the origin, z = 0, we 
set := (g*). 

Now, recall that the r-dynamics have the form 

= -V$(r) + 5(r)z 

where S is the 2 x N Stokes forcing matrix where the j-th column is given by the 
vector sin(kj • r) kj- /m | and {ki, . . . , kjv} is an enumeration of the active mode 
set K.. 

Since there is a positive probability that an Ornstein-Uhlenbeck process will 
remain in a small tube about any given continuous curve over some compact in- 
terval of time, as long as S remains non-degenerate (has rank 2), we may use the 
{zk} to drive the connector r to a neighborhood of r*. We shall later fix a 5^ < 6 
and denote 

n := mf{t > : |r(t) - r*| < 5r}. 

It may be that ||z(ti)|| is rather large. We must then bring the control to suffi- 
ciently fast so that at time 

T2 := inf{i > n : ||z(t)|| < 6^} 

we have |r(r2) — r*| < 26r- After this it must be shown that the process X can 
be held in place near the reference point x* through to the end of the designated 
interval /. The value for dz will be chosen to accomplish this last requirement. 

This argument may be extended to include general initial starting position due 
to the restriction to the compact region Q, where the coefficients of r are bounded 
uniformly. 

We now supply the details. 
Proof of Lemma \L4\ Construction of the norm-minimizing control: First, we dis- 
cuss the non-degeneracy of the Stokes matrix S. Note that some of the columns of 
S may become zero for certain r; however, by hypothesis, K, contains at least three 
pairwise linearly independent vectors, so S{r) will still have rank 2 for all r G Qy-. 

Let r(t) be a piecewise smooth curve in the interior of the annulus Qr with 
r(0) = ro and r(T) = r,,,. Without loss of generality, we will suppose that this 
path is a single linear segment since any two points in an annulus can be con- 
nected by two line segments. Furthermore, we suppose that T{t) is an arclength 
parametrization so that | J^r(t) | = 1. We will write T{t) = vt + tq where v is the 
unit vector ^(r,,, — tq). 
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For each fixed t G [0, T], the Unear system 

s{nt))z{t) = jT{t)+v^m) 

has a unique minimal norm solution, namely 

z(t) = 5+(vt + ro)(v + V«>(vt + ro)) (19) 

where 

= s*{ss*)-^ 

is the Moore-Penrose pseudoinverse IIBIG80II and S* is the transpose of S. Note 
that SS* is a symmetric 2x2 matrix and so its inverse can be written explicitly. 
Let Si denote the ith row of S. Then 

/QQ*\-l_f l'S'2p —Sl-S2 
[DD ) — I c- c 10 12 

The key observation is that each component is a sum of the form ^ cj sin^ (kj • 
(vi + ro)) where the coefficients cj G [—1, 1] do not depend on the path F. There- 
fore each component is a continuous function of t. This is true of S* as well and 
so z is in fact continuous. 

It remains to show uniform bounds on z, but this is immediate since V$ and 
S are continuous functions over the compact domain Qj.. Altogether we see that if 
N = \}C\,it follows from OH) that 

||z|| < sup |5"^(r)||l + V^>(r)| 

< |/C|2(l + max|V$(r)|). 

This is the constant seen in the definition of the good region, Eq. (fT/] ). 
Error bounds for fuzzy control: Let the r-path F be specified and let z be the 
associated norm-minimizing control. Let z{t) be a solution to the fluid velocity 
SDE (O and define to be the event where the sample path z(t; co) satisfies 

sup \z{t) — z{t; uj)\ < e . 

te[0,ri] 

We note that F{$7e} =: > 0. Now, let r{t;u!) be the solution to the uj-hy-u 
ODE, 

= -V$(f) + S(r)z . 

at 
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The control error h(i; io) := T{t) — r{t; uj) satisfies the ODE 

^|h|2 = -2(V$(r) - V$(f)), h) + 2(5(r)z - S{r)z, h) . (20) 

Since the spring potential <^ is continuously differentiable, it is Lipschitz in the 
region of good r-control Q,. with constant A$. The first term of (l20l ) is therefore 
bounded for every t by 

-2(V$(r) - V$(f)),h) < 2|V$(r) - V$(f))||h| < 2A$|h|2. 

The Stokes matrix is similarly Lipschitz with constant A5 and together we have for 
all u e fie 

^|h|2 < 2A$|h|2 + 2(5(r) (z - z), h) + 2{{S{T) - 5(f))z, h) 

<2|h|2(A$ + A5|z|) + 2e|5(r)||h| 

< |h|2(2A$ + 2A5|z| + e\S{T)\) + e|5(r)| 

where in the last line we have used the polarization inequality |h| < ^(1 + |hp). 
By Duhamel's principle we have 

|h(t)|2 < e rer(2A*+2Adz(s')|+e|5(r(s'))|)ds'|5(r(^))|^^ 
Jo 

We conclude that, given F, ri, and any 6r > 0, we may choose e,. sufficiently small 
so that 

pj sup \T{t) - r{t)\ <5r\>a,^> 0. (21) 

[t6[0,ri] J 

Settling the noise and holding r in place: Although f(Ti) is within a ball of 
arbitrarily small size about the target , it is likely that the magnitude of the forcing 
terms z will be nontrivial. We must show that it is possible to decrease the forcing 
terms rapidly enough that f does not leave a prescribed ball around r*, in this 
case radius 26r- Subsequently, we must hold r in place through to the end of and 
prescribed interval /. 

For alH G / with t > ri, we estimate as follows: 

\r{t) - r,| < \r{t) - f (rs)! + |r(r2) - ?(ri)| + |f (n) - r,| . 

where we recall that n is the arrival time of r in a neighborhood r* and T2 is the 
arrival of the full Markov process near x* . The estimates follow from the integral 
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representation of the f process. Over the interval [t2, t] we have 



|f(t)-f(T2)| = / -V^{^{s)) + S{r{s))i{s)ds 



<{t-T2) sup <^|V$(r)| + |5(r)| sup {||z(s)||} J 

reBs^ir,) [ se[r2,t] J 

Since the Stokes matrix satisfies the bound |S'(r)| < A^, we require ||z(s)|| < 
5/6Nt2 for all s G [T2,t\. Next note that since the spring potential $ is contin- 
uously differentiable and since r^, is located at a local minimum, we may choose 
6r G (0,(5/3) small enough so that sup{| V4'(r)| : r G S25,(r*)} < 

In the interval [ri , T2] we prescribe the ideal path for the fluid vector z to be the 
linear interpolation between the initial condition z(ti) and z(r2) = 0. We demand 
that the sample paths z{t) satisfy 

t — Ti 

\\z{t) — z{t)\\ < €r — {er — 



T2 - n 



for all t G (ti, T2) where tr is the constant chosen in the inequality (|2TI ) and T2 will 
be chosen in a moment. We have the estimate 

|f(T2) - f (ri)| < {t2 - Ti) sup [|V$(r)| + S{r) sup {||z(s)|| 
< (r2 - Ti) [ sup |V#(r)| + N{M + e)] 

where M is the constant in the definition ([TT] ) of the good set Q^, and the e allows 
for the tube surrounding the ideal path z. We may therefore choose a T2 sufficiently 
close to Ti so that |f (T2) — f (ti)| < 6/6. 

Altogether, we have for any r2 < t < i2, that |r(t) — r*| < |(5and ||z(t)|| < ^6 
with positive probability, so the claim ([TSl l is satisfied. 

□ 

We now use the properties of neighborhoods of the r singularity (and hence the 
"bad" set B) established in Corollary 1 1.3 1 to extend the previous proposition to all 
of C. 

Lemma 2.5 (Topological irreducibility on C). Given a 6 > 0, there exists at\> {) 
so that for any t > t'^ there is an a'^ > with 

xec 
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Proof. Set t\ = ti + 1 where ti is the constant from Lemma I2i4l Now for any 

t > 



inf Pt(x, Bsix,)) > inf V^{t < 1} inf inf P,(x, ^^(x,)) 
> aai > 

where a comes from Corollary 11.31 and qi from Lemma [241 Setting a'l = aai 
completes the proof. □ 



2.3 Measure Irreducibility via Hormander's Condition 

Lemma 2.6 (Absolute continuity of the transition density). For any t > 0, the 

Markov process {X(t) = (2(*))}t>0 with transition kernel 'Pt(x, U), possesses a 
transition density pt{x,y), i.e., 

Vt{^,U) = [ pt{x,y)dy 
Ju 

for every U £ ^{C), where pt{x,y) is jointly continuous in (x, y) £ C x C. 

Remark 2. In fact, the system has a density for all (x, y) G X x X. However, due 
to the periodicity of our forcing, proving this would require an additional small 
argument. Since we do not need this fact, we refrain. 

Proof. The claim follows from a now classical theorem of Hormander which states 
that if a diffusion on an open manifold satisfies a certain algebraic condition then 
Li = dt — C and L2 = dt — C* are both hypoelliptic in C where C is the generator 
of the diffusion X(t) and C* is its adjoint. A combination of Ito's formula and 
the fact that we have shown that the singularities of the potential are unattainable 
demonstrates that Liu = and L2U = have distribution- valued solutions. Hy- 
poellipticity of the operators ensures first that these distribution-valued solutions 
are in fact smooth. Furthermore, hypoellipticity implies the existence of funda- 
mental solution, which in turn yields continuity in the second variable throughout 
the center of the space C. 

The fact that the density is jointly continuous follows after a little more work. 
The argument is laid out in its entirety for valued diffusions in Section 7.4 
of llStrO Si In particular, see Theorem 7.4.3 and Theorem 7.4.20. Essentially, the 
same proofs follow in our setting since we have shown the system is a well defined 
diffusion on the manifold X with distribution-valued solution. Hypoellipticity and 
the properties which follow are local statements, and therefore still apply. The 
needed results in the general setting, as opposed to M^, can be found in Chapter 22 
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of HHorSSII . noting in particular Theorem 22.2.1. However, the presentation in 
llStrOSl is closer to the exact statements we need. 

We now turn to the explicit calculations needed to show that Hormander's con- 
dition is satisfied. We recast the system of equations ^ and ^ as 

dX{t) = A(X) + BdW{t) 
where A(x) G R^+^ and B G M(2+iV)x(2+7V) ^jj^ 

-V<I>(r) + EkSm(Ak.r)gzk\ „ ^ 0\ 

-x^uz J ' yo Ij ' 

Here / is the identity in M^. In this notation, the generator C of the diffusion is 
given in terms of a test function ip by 

(£^)(x) = {A • V)99(x) + IY.(^u- V)V(x) 

ke/c 

where is the column of B associated with the mode direction k G /C. Hormander's 
condition at the point x requires that 

spanjSk, [A(k), Sk] : k G /C} = R'^^^ 

where [^(x), i?k] is the commutator or Lie bracket between the two vector fields 
A(x) and i?k- In our simplified setting where i?k is a constant vector-field one has 

where ek is the is the unit basis vector in = M}^^ associated to the mode 
direction k G /C. 

The set { [^(x), i?k]}ke^ will span M2+^ if and only if the set {sin(k-r)k-L}ke/c 
spans since the {ek : k G /C} spans M^. We recall that by assumption /C con- 
tains at least three pairwise independent vectors which we label ki,k2, and ks. 
One may note that due to the periodicity of the forcing, sin(Ak • r) = for all 
r G 1? j}?. However, by construction, all of these points lie outside of C. Thus 
restricting to x G C at least two of r • kj are nonzero. □ 
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2.4 The Lyapunov function 

We now provide the Lyapunov function required by Theorem 12.21 to match the 
minorization established in Lemma [231 We only need the function to control the 
return to the center of space. As such we will truncate the function inside of a 
radius Rq. 

Recall the definition of the Lyapunov function 

y(x) := |r|2 Vi?g + r?||zf . 
where rj is specified below. 

Lemma 2.7 (Lyapunov function). There exists an interval Ii C (0, oo), a function 
y : X ^ [0, oo] with compact level sets, at least one point x such that V^(x) < oo, 
and a constant ci > such that for any t £ Ii there exists a cq £ (0, 1) with 



{rtV)i^)<coV{x) + ci. 

^r^acc + \ — ^^^^^^ 



Proof. Again, let C be the generator for the Markov process X(t) := (^[*|) , which 
we write explicitly as 



C := (- V<f (r) + U(r, z)) • Vr + ^X' ( E -^l^l'^k^ + 



Because the definition of V and estimates for the spring potential vary depend- 
ing on the value of |r| we divide the proof into three cases. First we address the 
case, |r| > Rq, where V$(r) • r > 7|rp. Denoting k := minkg/c{|k|}, we have 



£y(x) = -2V«>(r) • r + 2U(r, z) • r + rjuX^ ^ {-2\k\^zl + Pal) 

< -27|r|2 + 2 (^i?|U(r,z)|2 + i|r|2^ +r]uX^ {-2k^\\zf + f3\\a\\l) 

< -2 (7 - |r|2 + 2 {vi^X^^ - ^) \\zf + w>^'Ph\\l 

For 5 G (0, 1), letting = (1 — <^)7 and defining Ci := r//3z/A^||(T||Q yields 

£V{x) = -2(57|r|2 - {-quX^k'^ - [(1 - 6)-f]~^) ||zf + Ci . 

Further demanding that 6 < min(l, z/A;^A^/7) allows us to pick rj > [7(1 - 
5)(A2j/P-(^7)]-isothat 

£y(x) < -57y(x) + Ci 
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when |r| > Rq. 

When |r| < Rq, the r-dependent portion of the Lyapunov function is constant. 
Therefore 

CV{x.) < -2r]u\^k^\\'zf + Ci 

= -2v\^k'^{RQ + 7?||zf) + 2\^uk'^RQ + Ci 
= -2\^vk^V{^) + C2 

Let a = (57 A 2\^vk'^ then, taken together, the above two estimates yield for 
X(to) = xo with |r| / Rq, 

jE^, [V{Mm<-aV{Xt) + C2 

so that 

Exo[^(X(t))] < e-'^V(xo) + ^(1 - e-«*) 

Given t, we may set cq = e~"* and ci = 6*2/0 to give (fT2l) . 

It remains to include the case where X(to) = (^JJ) where |ro| = -Rq- We first 
observe that the measure irreducibility Lemma 12.61 guarantees that for any small 
h > 0, P{|X(to + h)\ = Rq} = 0. Furthermore, continuity of the dynamics and 
of the Lyapunov function ensures that there exists an e > such that if we define 
X(t) = (f(t),z(t)) with r{to) = {Rq + e)X(to)/^o and z(to) = z{to), then we 
have 

VCKito + s)) < V{X{to + s)) 

for all s < h. 
It follows that 

Exo[l^(X(t))] < e-'^V(xo) + ^(1 - e--'). 
and taking limits as h and e go to zero, we recover ([T2l ). □ 

2.5 Ergodicity of generalizations 

In the derivation of the model equations ([5]) and® we imposed the simplifying 
assumption that the center of mass m(t) := i(xi(t) + X2(t)) is held at zero (see 
Appendix). This greatly simplified the presentation, and did not affect the conclu- 
sion that the bead-spring system has an ergodic connector process r{t). Indeed the 
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fluid velocity term with nonzero m(t) is given by Eq. (|27] |: 

^[U(xi(t),t)-U(x2(t),t)] 

= [cos(Ak • m)2;k — sin(Ak • m)yk] sin(Ak • r)-|j^ 

where the {yk} are a second set of OU-processes defined exactly as the {^k}- 

Because the m terms appear inside of cosines and sines, there is no new signif- 
icant contribution to the Lyapunov function calculation. For the Hormander con- 
dition, the additional terms in the coefficients of the noise introduces more "dead 
spots" in the forcing, but still one needs only four pairwise linearly independent 
vectors kj in the mode set JC to ensure that at least two of the vectors 

{ [cos(Akj • m) — sin(Akj • m)] sin(Aki • r)k^} 

are nonzero. This guarantees the existence of a continuous transition density and it 
remains to show the (5-ball controllability as in Lemma l2!4l While the calculation is 
more involved, the principle of identifying the region of good control Q, where the 
coefficients of the r-differential equation are uniform, still applies. Furthermore, 
since the differential equation for r is linear in the {yk} and {z]^}, we may still 
solve for stochastic control explicitly in terms of the desired path T as long as the 
new Stokes matrix is non-degenerate. Again, this is guaranteed by the hypothesis 
that JC contains at least four pairwise linearly independent vectors. 

We take a moment to consider the model closest to that of Celani, et. al. 
IICMV05II . where in the canonical Langevin Equation ([U, the spring potential is 
quadratic, the mass m is still 0, but the coefficient of the Brownian motion is 
nonzero: k = -v/SfceT^. Our generalization is the replacement of the Kraichnan- 
ensemble with the stochastic Stokes equations. 

Again, for small values of r the force on the connector r due to the fluid veloc- 
ity becomes negligible, however, the remaining terms an Ornstein-Uhlenbeck and 
by standard ergodic properties of such processes, r quickly leaves any neighbor- 
hood of the origin with probability 1. For large values of |r|, the quadratic spring 
potential dominates and the Lyapunov function calculation still holds. Since the 
diffusion is elliptic, trivially implying the existence of a continuous transition den- 
sity, and all arguments in the derivation of the stochastic 5-ball controllability still 
apply, we see that the ergodic theorem holds for r(t). 

This stands in contrast to the results in IICMV05II where it was argued that 
there exists a range of parameters where no stationary distribution exists. It is not 
immediately clear to us how to construct a model with coloured-in-space-and-time 
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fluid velocity field that supports the "stretched" and "coiled" regimes cited in the 
physics literature. 
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A Derivation of the model 

In the overdamped, highly viscous regime, it is reasonable to neglect the nonlinear 
term in Navier-Stokes equations ||OR89i Following ||Wal86L IIDZ92L ||Dal99ll and 
IIMcK06ll we have the stochastic PDF given in Section 1, Fq. |2l 



with periodic boundary conditions on the rectangle R := [0,L] x [0, L] where L 
is presumed to be very large. We assume that the space-time forcing is a Gaussian 
process with covariance 



where a, /3 € {xi, X2} and 6a/3 is a Kronecker delta function. The constants each 
have physical meaning: ks is Boltzmann's constant, T is the temperature of the 
system and u is the viscosity. It follows that 



where the coefficients a"k are related to the spatial correlation function T through 
the Fourier relation 



(9iu(x, t) - z/Au(x, t) + Vp(x, t) = F(x, dt) 
V • u(x, dt) = 



E F°(x, t)F^{y, s) ={tA s)2kBTu6^pT{^ - y) 




kGZ2\0 




fcez2\{o} 
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E 



Indeed, from the definition ^ we have 

k j 



demonstrating ([3]l- 

We turn our attention to the Fourier transform of the SPDE, noting the trans- 
form of the noise is given by 



g— Aik-x 



jez2\o 

= y/2kBTiyL Y ^jBj(t)5kj 

jez2\o 

= v^2fc^LcJkBk(t) 
The SPDE transforms into the infinite dimensional system 



R 



duk(t) + A^z/|kpUk(i) + Azkpk(t) = y^2k^fuLakdBk{t) (22) 

Aik • u(t) = (23) 

For the sake of completing the formal argument, suppose for the moment that the 
forcing term is smooth with derivative f . By taking the dot product of k with the 
terms of equation (l22l ). the first two terms vanish - via incompressibility condition 
(|23] ) - leaving the identity 



Ai|k|2pk(t) = \/2fc^IVLcrkk • i{t) (24) 
Substituting back into (l22l) and gathering f (t) terms on the right-hand side yields 

duk(t) + A2zy|k|2uk(t) = y^2kBTuLai, (^f (t) - ^J^^^ (25) 
The projection on the right hand side has two standard representations. 



24 



where k-*- := (^^^) . Applying Duhamel's principle and assuming vanishing initial 
conditions, we have the following representation for solutions to the fluid mode 
equations 



Zk 



|k|2 

where we define to be the appropriate Omstein-Uhlenbeck process, 

dzk(f) = -X^u\kfz]^{t)dt + \/2kBTuLakdBk{t) . 
We therefore have the solution for the fluid velocity field, 

kez2\o ^ 11/ 

L2 - |k|2 

kez2\o ' ' 

After defining ■= ' have the l-d OU-processes that drive the dy- 

namics 

dz]^[t) = ^2 Zk[t)dt H dBy^{t) 

Imposing the condition that we require real-valued solutions, after Fourier inver- 
sion we have the following trigonometric expansion for 2-d stochastic Stokes. 

^ k-L 
u(x, t) = 2^ (cos(A k • x)yk + sin(A k • x)zk) -r— (26) 

kez2\o ' ' 



where the t/k and Z]^ satisfy 



dy.it) = -'-^y.it)dt + ^^^^f^dwyn) 



, , , 47r^j^|kP , , , -v/2^BTV(7k ^ 
dzkit) = jY^Zk{i)dt + -dW^it) 

and {W^ki {^kl are i.i.d. sequences of standard l-d Brownian motions. In 
the bulk of the paper we will express the above SDEs in terms of the constants 
A = 27r/L and /3 = ksT /2tt. 
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In the case where = 1 for all k, we have spatially white noise forcing. How- 
ever, in dimensions two and higher ||Dal99ll l|Wal86ll . there do not exist function- 

2 

valued solutions for stochastic Stokes. It is sufficient to assume < the 

condition specified in Assumption [T] (Note that the sum does not converge with 
(Tfc„ = 1 because J2]s.ez'^\o '^^e equivalent of a double integral.) 

In this paper, we study the dynamics of the two beads in normal coordinates: 

m{t) = i(xi(t) + X2(t)) and m(t) = i(xi(t) - X2(t)), 

^m(t) = i[u(xi(t),t)+u(x2(t),t)] 

|r(i) = -V$(|r(t)|) + ^[u(xi(t),t) - u(x2(t),t)]. 

In light of equation Q, we may write the radial process and the noise together as 
a Markovian system of SDE with two degenerate directions. In order to write the 
system in this form, we first record the identity 

^[u(xi(t),t) -u(x2(t),t)] (27) 
= [cos(Ak • m)zk — sin(Ak • m)yk] sin(Ak • t^)t^ 

For the majority of the paper, we used the simplification m(t) = for all t. This 
does not have any effect on the ergodic results as is discussed in Section 1231 but it 
does significantly streamline the presentation. Altogether we have the definition of 
the dynamics given in Section 1 , Eq. [5] 
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